DOI:  10.15415/jnp.2013. 11003 


Cluster  Radioactivity  in  127I 

M.  Balasubramaniam1  and  K.  Manimaran2 


'Department  of  Physics,  Bharathiar  University,  Coimbatore-641046,  India 

department  of  Science  and  Humanities,  Theni  Kammavar  Sangam  College  of  Technology, 

Theni-625534 

Abstract  Using  the  preformation  cluster  model  of  Gupta  and  collaborators  we  have  studied  all 
the  possible  cluster  decay  modes  of  127  I.  The  calculated  half-lives  are  compared  with  recently 
measured  lower  limits  of  cluster  decay  half-lives  (for  the  clusters  like  24Ne,  2sMg,  30Mg,  32Si,  34Si, 
48Ca  and  49Sc)  of  127I.  Our  calculated  half-life  values  lies  well  above  the  experimentally  measured 
lower  limits  and  the  trend  of  the  values  also  matches  with  experimental  ones. 
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1.  INTRODUCTION 


A  new  kind  of  exotic  radioactive  decay  mode,  in  between  a- decay  and  spontaneous 
fission  called  cluster  decay  or  cluster  radioactivity  or  heavy  particle  decay  apart 
from  the  three  basic  decay  modes  (a-decay,  /3-decay  and  7-emission)  was  predicted 
theoretically  in  1980  [1].  Later  in  1984  Rose  and  Jones  [2]  observed  experimentally  the 
emission  of  14C  from  223Ra  with  branching  ratio  of  (8.5±2.5)xl0-10  relative  a-particle. 
Subsequently  other  authors  [3,  4]  also  observed  the  same  cluster  14C  from  same 
radioactive  parent  223Ra.  A  few  years  after  the  above  observations,  several  other  decay 
modes  like  200, 23F,  24,26Ne,  28-30Mg  and  32,34Si  from  different  radioactive  parent  nuclei 
like  221Fr  22 1-224,226^  225 228,230,232^  231Pa  230,232-236^7  237jqp  236,238,240pu  241  an(j 

242Cm  were  observed  by  various  experimental  groups  in  all  around  the  world  with 
branching  ratio  relative  to  a-decay  from  10-9  to  10-16  [5,  6],  All  the  emitted  clusters 
are  heavier  than  a-particle  but  lighter  than  lightest  fission  fragment  observed.  The 
daughter  nuclei  observed  are  always  double  magic  nucleus  208Pb  (N=126  and  Z=82) 
or  its  neighbouring  nuclei,  which  implies  that  cluster  decay  process  associated  itself 
with  the  closed  shell  behaviour  of  emitted  daughter  nucleus.  Simultaneously  it  has  also 
been  studied  extensively  using  various  models  after  its  experimental  verification. 


In  general  there  exists  two  kinds  of  models  for  explaining  the  observed  decay  modes 
and  for  predicting  new  decay  modes.  In  one  kind  of  model,  the  a-particle  as  well  as  the 
heavy  cluster  (or  clusters)  was  assumed  to  be  pre-born  in  a  parent  nucleus  before  they 
could  penetrate  the  barrier  with  the  available  Q-value.  These  models  are  in  general 
called  as  ’’Preformed  Cluster  Models”  (PCM)  [7-11].  In  such  a  model,  clusters  of 
different  sizes  are  considered  to  be  preformed  in  the  parent  nucleus  with  different 
probabilities.  In  the  other  model,  only  Gamow’s  idea  of  barrier  penetration  is  used  i.e. 
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Balasubramaniam,  M.  without  considering  the  cluster  being  or  not  being  preformed  in  the  parent  nucleus. 

Manimaran,  K.  These  models  are  in  general  called  as  ’’Unified  Fission  Models”  (UFM)  [12-16]. 

Several  semi-empirical  formulae  were  also  proposed  to  calculate  the  partial  half- 
lives  of  cluster  decay  modes  in  trans-lead  region.  A  semi-empirical  formula  with 
only  three  parameters  to  calculate  logarithm  of  half-lives  of  cluster  decay  modes  was 
proposed  by  one  of  us  [17].  Recently  a  scaling  law  has  been  given  by  Horoi  et  al. 
to  calculate  the  logarithm  half-lives  cluster  decay  modes  [18].  Apart  from  the  trans¬ 
actinide  region,  it  was  theoretically  predicted  that  the  trans-tin  region  as  fertile  region 

26  to  observe  the  heavy  particle  decay  due  to  the  closed  shell  behaviour  of  Sn  nucleus 

[19-25].  Based  on  different  theoretical  models  Ba  and  Ce  isotopes  were  predicted 
as  cluster  emitters  to  emit  12C  and  160  clusters  respectively  leaving  Sn  as  daughter 
nuclei.  Later  12C  emission  is  reported  experimentally  from  112Ba  nucleus  and  from 
114Ba  with  upper  limit  for  the  half-lives  as  >  3.63  s,  1.70xl04  s  and  >4.10  s.  Recently 
a  new  semi-empirical  formula  is  proposed  by  us  [26]  which  is  a  modified  form  of 
[17]  to  calculate  the  logarithm  half-lives  of  12C,  160, 20Ne,  24Mg  and  28Si  clusters  from 
various  isotopes  of  Ba,  Ce,  Nd,  Sm  and  Gd  respectively.  Very  recently  Bernabei  et  al. 
[27]  measured  the  new  upper  limits  of  the  half-lives  of  24Ne,  28Mg,  30Mg,  32Si,  34Si, 
48Ca  and  49Sc  cluster  radioactivity  in  127I  parent  nucleus.  In  the  present  work  we  have 
theoretically  studied  all  the  possible  cluster  decay  modes  of  127I  using  PCM  proposed 
by  R.  K.  Gupta  and  co-workers  [7-9]. 


2.  PREFORMATION  CLUSTER  MODEL 

Preformation  Cluster  Model  (PCM)  of  R.  K.  Gupta  and  co-workers  [7-9]  is  developed 
by  adapting  the  Gamow’s  theory  of  o-decay  but,  instead  of  square  well  potential,  a 
more  realistic  nuclear  interaction  potential  is  used  and  also  a  preformation  probability 
Pfl  is  associated  with  the  size  of  the  cluster.  The  clusters  of  different  sizes  are  having 
different  preformation  probabilities,  which  decreases  with  the  increasing  size  of 
cluster.  Thus  the  half-life  and  decay  constant  (X)  in  PCM  is  defined  as 

T —  and  \  =  P0P.  (!) 

X  A 

with  P  is  the  barrier  penetration  probability,  PQ  is  the  preformation  probability  and  v  is 
the  assault  frequency  calculated  as  in  [7-9].  In  PCM  of  R.  K.  Gupta  and  co-workers 
the  preformation  probability  is  a  theoretically  calculated  quantity  by  solving  the 
Schrdinger  equation  of  the  motion  in  mass  asymmetry  (charge  asymmetry)  coordinate, 
at  fixed  R  (defined  later)  and  is  essentially  based  on  the  nuclear  structure  information 
of  the  decay  process.  The  mass  and  charge  asymmetry  is  defined  as 

A  —  A  Z.-Z, 

T)  = -  and  ?7  = -  (2) 

A  z  Z 
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where  A  and  Z  are  the  mass  and  charge  numbers  of  the  parent  nucleus.  A  and  Z  with  Cluster 

i=  1 ,2  corresponds  to  the  mass  and  charge  numbers  of  daughter  and  cluster  respectively.  Radioactivity  in  127I 
The  preformation  probability  P(|  is  obtained  by  solving  the  following  Schrodinger 
equation 


h1  d  1  d 
2^[b~  dt]  .[b~  dr) 

v  rm  v  w 


+  VR(v) 


Ip"  (v)  =  E"i)"  (jy), 


(3) 


Using  a  similar  equation  such  as  shown  above  one  can  solve  the  penetration  probability 
P  in  R  co-ordinate  for  fixed  r)  as 


h2  d  l  d 

2  Jb~  dR  .[b~  dR 

V  RR  V  RR 


+  Vr(R) 


rp"(R) 


E"rp"(R), 


(4) 


27 


with  R=C|+Cr  C.(i=l,2)  is  the  Sussman’s  central  radius  which  is  related  to  the  effective 
sharp  radius  R.  by  (b~  1  fm) 


C  =  R 


.-4' 

R 


(5) 


where  the  sharp  effective  radius  is  given  as  R  =  1.28A2/3-0.76+0.8A“1/3.  The  collective 
potential  or  the  fragmentation  potential  V(r|)  appearing  in  Eq.3  is  calculated  as  at  fixed 
eta  in  R-coordinate  (R-motion). 


y  (77)  =  Bt  (A ,  Z. )  +  yc  +  vp  (6) 


where  B(i=l,2)  corresponds  to  the  binding  energy  of  the  daughter  and  cluster 
respectively.  Vc  and  Vp  are  the  coulomb  and  proximity  potential  calculated  as  in  [7-9]. 
The  scattering  potential  appearing  in  Eq.4  is  simply  the  sum  of  coulomb  and  proximity 
potential  calculated  at  fixed  eta  in  R-coordinate  (R-motion) 

The  Q-values  of  the  decay  modes  are. 


Q  =  M  —  mi 


(7) 


where  M  is  the  mass  excess  of  the  parent  nucleus  and  m.  with  i  =  1 ,2  corresponds  to  the 
mass  excesses  (taken  from  [29])  daughter  and  cluster  respectively,  expressed  in  MeV. 
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Figure  1:  The  calculated  fragmentation  potential  using  Eq.  6  for  127I.  The 
clusters  with  local  minima  are  labelled. 


3.  CALCULATIONS  AND  RESULTS 

In  the  present  work  we  have  studied  the  all  possible  cluster  decay  modes  in  127I. 
Figure  1  presents  the  charge  minimized  fragmentation  potential  from  Eq.  6 
According  to  quantum  mechanical  fragmentation  theory  the  minimum  in  the 
fragmentation  potential  gives  rise  to  the  maximum  in  the  probability  to  observe  the 
particular  decay  mode.  The  fragmentation  potential  increases  with  the  increasing 
mass  number  of  the  clusters,  but  there  are  some  local  minima  in  the  potential  energy 
surface  for  the  clusters  like  4He,  10Be,  14C,  20O,  24Ne,  28Mg,  48Ca,  and  49Sc  in  which 
10Be,  14C,  200, 24Ne,  28Mg  clusters  are  already  observed  in  trans-lead  region.  Also  there 
exist  a  small  cold  valley  in  the  near  symmertic  region  and  having  minimum  in  the 
potential  double  magic  cluster  48Ca  (N=28,  Z=20),  and  its  neighbour  nucleus  cluster 
48Ca  (N=28,  Z=20),  and  its  neighbour  nucleus  49Sc.  Figure  2  presents  the  calculated 
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Cluster 
Radioactivity  in  127I 


Figure  2:  The  calculated  Q- values  of  cluster  decay  modes  of  127I  using  Eq. 
7.  The  Q- values  are  positive  only  beyond  the  cluster  with  mass  number  22. 
The  clusters  with  local  maxima  are  labelled. 


Q-values  which  is  the  available  energy  for  the  clusters  to  penetrate  the  potential 
barrier.  The  maximum  in  the  Q-value  increases  the  penetration  probability  which 
increases  the  probability  to  observe  the  particular  decay  mode.  Though  the  smaller 
clusters  like  10Be,  14C,  20O,  24Ne  has  minimum  in  potential  energy  their  Q-value  are 
very  low.  The  Q-value  systematics  prefers  the  heavier  clusters  (48Ca,  and  49Sc)  in  near 
symmetric  region  to  observe  such  cluster  radioactivity. 

The  Penetration  probability  is  calculated  only  for  the  clusters  having  positive 
Q-values  and  is  presented  in  Figure  3.  The  clusters  like  26Mg,  43K,  46Ca,  50Ti,  53V, 
and  5SCr  has  the  maximum  in  the  penetration  probability.  Figure  4(a)  presents  the 
calculated  preformation  probability  P0  for  the  complete  mass  asymmetry  involved, 
as  mentioned  earlier  the  preformation  probability  decreases  with  increase  of  mass 
number  of  the  clusters.  There  exist  large  fluctuations  and  increase  in  the  PQ  for 
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Figure  3:  The  calculated  penetration  probability  (P)  only  for  the  cluster 
decay  modes  l27I.  having  positive  Q-values.  The  clusters  having  minima  in 
P  are  labelled 


the  clusters  in  the  near  symmetric  mass  region  and  has  maximum  probability 
for  the  48Ca  (N=28,  Z=20),  and  its  neighbour  nucleus  49Sc.  For  clear  vision  of 
the  fluctuations  near  symmetric  mass  region,  it  is  enlarged  and  presented  in 
Figure  4  (b). 

The  decay  constant  and  logarithm  of  half-lives  all  possible  cluster  decay  modes  of  127I, 
with  positive  Q-values  are  calculated,  presented  in  Figure  5  and  6  respectively.  It  is 
clear  from  the  figure  5  the  decay  constant  reflects  the  combined  effect  of  penetration 
probability  and  preformation  probability  since  the  assault  frequency  is  merely  a 
constant  varies  in  between  102°  and  1022.  The  decay  constant  prefer  the  clusters  that 
posses  maximum  in  both  the  penetration  probability  and  preformation  probability.  The 
log  T1/2  values  also  prefer  the  same  clusters  that  are  preferred  by  the  decay  constant 
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Figure  4:  (a)  The  calculated  preformation  probability  (P0)  for  the  whole 
mass  asymmetry  involved  in  127I.  (b)  The  rectan-gularly  marked  portion 
in  (a)  is  enlarged  for  clear  vision  of  preformation  probability  (P0)  for 
clusters  in  near  symmetric  region.  The  clusters  with  maximum  PQ  and  their 
complementary  daughters  in  this  region  are  labelled. 


31 


since  log  Tlfl  is  derived  from  the  decay  constant.  The  clusters  (24Si,  28Mg,  30Mg,  32Si, 
34Si,  48Ca,  and  49Sc)  for  which  limits  of  cluster  decay  half-lives  is  measured  in  the 
recent  experiment  [27]  is  labelled  in  both  figure  5  and  6.  Figure  7  presents  comparison 
of  our  calculated  cluster  decay  half-lives  with  the  measured  lower  limits  of  cluster 
decay  half-lives  of  127I.  Our  calculated  half-lives  lies  well  above  the  experimental 
lower  limit  for  the  clusters  30Mg,  32Si,  34Si,  48Ca,  and  49Sc.  The  trend  of  the  variation  of 
cluster  decay  half-lives  for  heavier  clusters  (32Si,  34Si,  48Ca,  and  49Sc)  matches  with  the 
measured  half-lives. The  calculated  half-lives  lies  below  the  measured  lower  limit  for 
clusters  24Ne  and  28Mg  and  this  may  be  due  to  the  very  lesser  Q-values. 
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Figure  5:  The  calculated  decay  constant  (X)  using  Eq.  1  for  the  possible 
cluster  decay  modes  of  127I,  with  positive  Q-values.  The  clusters  for  which 
the  experimental  lower  limits  of  half-lives  measured  in  [27]  are  labelled. 


Fragment  mass  number  (A2) 

Figure  6:  The  calculated  log  values  of  half-life  for  the  possi-ble  cluster 
decay  modes  of  127I,  with  positive  Q-values.  The  clusters  for  which  the 
experimental  lower  limits  of  half-lives  measured  in  [27]  are  labelled. 
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Figure  7:  The  calculated  logarithm  of  half-lives,  expressed  in  seconds,  are 
compared  with  the  experimentally  measured  [27]  lower  limits  of  possible 
cluster  decay  modes  of  127I. 
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Table  1:  The  table  shows  log  values  of  our  calculated  cluster  decay  half-lives  using  PCM  and 
log  values  of  experimentally  measured  lower  limits  of  cluster  decay  half-lives  of 127 1. 


Cluster  decay  modes 

'°glOTl/2(S) 

present 

Ref.  [27] 

I27I->  24Ne  +  103Tc 

20.07 

30.65 

127I^  28Mg  +  "Nb 

26.10 

29.80 

127I — ^  30Mg  +  97Nb 

31.58 

31.82 

127I^  32Si  +  95Y 

33.23 

28.98 

127I^  34Si  +  93Y 

36.10 

30.24 

127I^  48Ca  +  79 As 

33.28 

29.33 

127I^  49Sc  +  78Ge 

33.18 

28.95 
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4.  SUMMARY 

Summarizing  in  this  work,  we  have  studied  the  cluster  decay  half-lives  of  all  possible 
cluster  decay  modes  of  127I  nucleus  using  preformation  cluster  model.  The  calculated 
cluster  decay  half-lives  of  different  clusters  (24Si,  28Mg,  30Mg,  32Si,  34Si,  48Ca,  and  49Sc) 
emitted  from  127I  are  compared  with  the  experimentally  measured  lower  limits.  Our 
calculated  cluster  decay  half-lives  lies  well  above  the  experimental  lower  limit  except 
for  two  lighter  clusters  (24Si  and  28Mg).  The  trend  of  the  calculated  values  also  matches 
with  the  experimental  values  for  the  heavy  clusters  30Mg,  32Si,  34Si,  48Ca,  and  49Sc. 
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